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, , Abstract 

X 

5^ I We study QCD hard processes at scales of order ^ in the limit in 

which the beta-function coefficient b is taken to be small, but as(k^) is kept 
fixed. The (nonperturbative) Pomeron is exponentially suppressed in this 
limit, making it possible to define purely perturbative high-energy Green's 
functions. The hard Pomeron exponent acquires diffusion and running cou- 
pling corrections which can be expanded in the b parameter and turn out 
to be dependent on the effective coupling ba'^Y. We provide a general setup 
for this 6-expansion and we calculate the first few terms both analytically 
and numerically. 



1 Introduction 



High energy hard scattering has received considerable attention in recent years, due to 
the increase with the squared center-of-mass energy s of the experimental cross section 
seen at HERA The analysis of high energy QCD shows that, at leading logs 
level in which as is considered to be frozen, the cross section is power behaved, and its 
exponent cug = asXm = ^^Xm is provided by the saddle-point Xm of the characteristic 
function x{l) of BFKL kernel. 

When higher order corrections are taken into account 0, the theoretical analysis 
changes conceptually, because of running coupling effects and of higher loop contribu- 
tions to the kernel. First of all, 0;^ = as{t)xm becomes t-dependent (where t = Infc^/A^ 
and k is the momentum of the hard probe) and is no longer related to the leading sin- 
gularity in the cu-plane. Furthermore, the fact that as{t) becomes larger at small values 
of k"^ = (9(A^) causes the existence of bound states (isolated, or in some cases dense 
u singularities) so that the leading one, the Pomeron tup, becomes crucially dependent 
on the strong coupling region and thus on the soft physics. The question then arises of 
what one can really compute and measure in a perturbative QCD approach. 

It is a common belief that the "hard Pomeron" exponent is actually measurable in 
two-scale hard processes (for example the two jet production in hadronic collisions [Q, 
the forward jet/vr^ production in deep inelastic scattering p or the 7*(A;)7*(A;o) X 
process 0) in which for sufficiently large A;^, fcg ^ A^, the cross section is roughly 
determined by the exponent a;^, evaluated at some average scale kko. It is also believed, 
though, that diffusion corrections ^, ^ to this exponent are to be expected, and that 
for sufficiently large Y = In ^ the asymptotic, nonperturbative Pomeron tup takes 
over. 

The interplay of the hard Pomeron regime (large fc^'s) with the Pomeron one (large 
Y) has been the subject of intensive studies p!0| , |1T[, mostly dealing with running 



coupling effects and diffusion corrections. Recently, it has been pointed out [0, |T^ for 
various BFKL-type models that the transition to the Pomeron regime is more like a 
sudden tunneling phenomenon rather than a slow diffusion process. A key result of this 
analysis is that for /c ^ /cq ^ A, the hard Pomeron behaviour dominates if uJs{kko)Y is 
large but not too large, its critical value being of the order of i, if uj-p ~ ^ parametrises 
the size of the Pomeron exponent {t is the scale at which ds{t) is regularized). Thus, 
the puzzling question arises again: is there really room for an observable hard Pomeron 
regime? 

Some source of optimism comes from the fact that the physical Pomeron of soft 
physics is indeed rather weak (the exponent is around 0.08 in most estimates fl^). 



Thus if LUjp is such a Pomeron, the effective critical value of UsY may be sufficiently 
large for the hard Pomeron to be seen. It is true that, the physical Pomeron is a 
complicated unitarity effect which turns out to be weak and approximately factorized 
for well studied but not well understood reasons. However, the trend of subleading 



corrections ^ is that of reducing the high-energy exponents and the diffusion coefficient, 
and this presumably slows down the onset of the nonperturbative behaviour as well. 
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A different way of studying the problem, motivated in the present paper, rehes 
on a somewhat different, but crucial observation (Sec. 2). The size of the tunneling 
probability at scale t = Infc^/A^ is easily seen to be of the order exp[— ^^^(/((^^(t))], 
where b is the beta-function coefficient which determines the running of the coupling 
as{t) = l/bt, and g{as(t)) is a function whose exact form depends on the nature of 
the problem being studied. Therefore, in the formal limit 6 — with taken to be 
fixed there is no tunneling to the Pomeron and only the hard Pomeron survives. In this 
paper, we take advantage of the Pomeron suppression for 6 ^ in order to study the 
perturbative part of the gluon Green function by expanding it in the b parameter both 
analytically and numerically. It turns out, that in BFKL-type models b plays the role 
of /i in a semiclassical expansion, and thus the above observation is actually the basis 
for a full calculation of the Green's function exponent which allows to find the diffusion 
corrections to the hard Pomeron in a systematic way (Sec. 3). 

We stress the point that we are expanding in b the gluon Green function with 
two hard scales {t,to ^ 1) by keeping as(t) and possibly as(to) fixed, we are not 
expanding the renormalization group logarithms themselves. Therefore, our expansion 
is somewhat different from previous related expansions of the Green's function ^, ^ 
and from analogous expansions applied to the collinear factorization limit [|I^, [T^ . We 



also stress that our expansion concerns the perturbative part of the gluon Green's 
function, the Pomeron one being exponentially suppressed as stated before. Here a 
technical point arises: how to define such perturbative part? We argue in Sec. 3 that 



the customary [10| double 7 representation can be given a precise perturbative meaning 



in a general BFKL-type model, apart from the intrinsic renormalon ambiguities, which 
come from the Landau-pole regularisation and are suppressed like ~ exp(— ^^i^). We 
thus overcome the difficulties found in Ref. [|] for the principal value prescription, by 
a careful choice of the integration contours. 

The 6-expansion provides a hierarchy of diffusion corrections to the Green's function 
exponent for t ~ to depending on the parameter z = |6^a^XmXm'^s^^; introduced in 
Ref. [|l^, and which turns out to be of the order 6^. For z <^ 1 or UsY <C t we recover 



the well known 0, |], |[ a^Y^ terms, and we compute new ~ af^^ and a^Y diffusion 
corrections which change the normalisation and the exponent at order 6^ (Sec. 3). The 
estimate of the Y^ and Y'^ terms is confirmed by the numerical methods (Sec. 5), which 



could be extended to the realistic model [IS] with the purpose of resumming sub leading 



effects, as in the improved equation or in the alternative approaches ||2T], 



41 25 1 . We are also able to evaluate the convergence radius Zc of the above series. 



and the large z behaviour which turns out to be oscillatory |jT8| and is actually masked 



by the onset of the asymptotic Pomeron. An amusing result is an essentially explicit 
evaluation of the perturbative Green's function for the Airy diffusion model (Sec. 4) 
which improves the estimates of Ref. ||18|| . 

Another application of the 6-expansion is in the numerical determination of high- 
energy Green's functions. Current practice in numerical solutions of the small-x evo- 
lution equations (see for example [2^, ^) is to evaluate a perturbative result with an 
infrared regularised coupling (e.g. with a cutoff) and then define the perturbatively 
accessible region in Y as that insensitive to the choice of regularisation. While physi- 
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cally intuitive, this approach has two significant disadvantages. Firstly the choice of a 
particular set of regularisation schemes is quite subjective. Secondly the formal limit 
on the perturbative region is determined by tunneling, i.e. F^ax ~ ^/o-sito)^ whereas 
perturbative predictions can in principle be made up to values of Y of order l/a^(to)- 
As is discussed in section the 6-expansion offers a solution to both these problems. 

The overall picture emerging from this work is that the hard Pomeron and its 
corrections can be properly defined in the 6-expansion, and its features can be found 
either numerically or analytically at semiclassical level in a precise way. Beyond this 
level of accuracy some ill-defined exponentially suppressed terms come in which are 
nonperturbative in both h and cxs and among them the asymptotic Pomeron itself. It 
is not known whether the size of such contributions may be hinted at or restricted on 
the basis of perturbative arguments. 

2 Hard Pomeron vs Pomeron in the perturbative 
regime 

The basic problem of small x QCD is to determine the gluon Green's function G{Y; t, to), 
where the logarithmic variables Y = log s/kk^, t = log k^A^ and = log ko^/A^ are 
related to the center-of-mass energy ^/s and to the transverse momenta of the gluons k, 
ko {k = |k|, kfj = |ko|) and A is a QCD parameter. In BFKL-type models, the function 
G satisfies the following evolution equation 

^G{Y;t,to) = K®G 

G(0;t,to) = 5(t-to), (1) 

where is a kernel in t-space. It is well known p| that, at leading level in logs, the 
kernel K{t,t') = asKo{t — t'), is scale invariant with symmetric characteristic function 
(^sXoil) = «sXo(l — 7), so that the high-energy exponent of Eq. (|l]) becomes 

ujjp = uj°, = a^max xo(^ + = «sXo(^), = (2) 

o-eR / Z TV 

Therefore, in the LLA this exponent is both the leading singularity in the u plane 
(conjugated to Y) and the saddle point value of the characteristic function x- For not 
too large values oft — to, a good approximation of G is provided by the diffusion model 

G{Y;t,to) ^ , ^ exp ( u;',Y - ^-^^^) , = l^iilM . (3) 

When subleading corrections in logs are evaluated running coupling and higher 
loop corrections come into play. The kernel K is no longer scale invariant — being 
dependent on a sit) in a essential way — and may have, besides the continuum spectrum, 
a discrete one whose maximum eigenvalue is uj-p. Furthermore, the fact that as{t) 
becomes stronger at small values of k'^ ~ 0{h?) favours the diffusion towards small 
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t values, so that eventually the asymptotic exponent in Y, that is tup, is crucially 
dependent on the soft physics. The magnitude of tup will then depend on the details of 
the physical regularisation of dsit) adopted in the region around the Landau pole. We 
shall refer for definiteness to kernels of the slightly asymmetrical form|^ 

K = a,{t)K{t-t') (4) 

where as{t) embodies the strong coupling boundary condition below some scale i and 
K is scale invariant. 

On the other hand it is still believed that, when t and to are fixed by the experimental 
probes to be sufficiently large (for example in scattering of two highly virtual photons) 
the diffusion to low t's is suppressed and the s dependence is basically power-like and 
t-dependent 

G{Y;t,to) ~ , ^ ^=exp + corrections ) , Y> t,to >1, 



(5) 

where the exponent uJs{t) is now perturbatively calculable. It should be noted however 
that uJs{t) is no longer a singularity in the u plane in that case. 



Some insight into the interplay of uJs{t) and ujjp is gained [|T8[ in simplified models 



like the Airy diffusion model or collinear models |jT2| in which the gluon Green's 



function in u space is factorized for a kernel K of the form (^) as follows 

GUt,to) = ^-^ {fFUt):F^ito)Qit - to) + t< — .to) , (6) 

where J^uj{t) is the (right) regular solution for t — > cx) of the stationary homogeneous 
diffusion equation {u — K)T = 0, and to-^(^o) is instead the (left) regular one for 
to — i> —CO {ko — * 0). In such models the strong coupling boundary condition on 
determines it in the form 

^-(^o) = ^'(to) + S{u)TM, (7) 

where J-'^ is now an irregular solution for to oo and the reflection coefficient >S'(ti;) 
carries the leading spectral singularities of the problem, in particular the Pomeron 
singularity wp. By using the decomposition (^ of JF into its perturbative and nonper- 
turbative components in representation (|^) we obtain the following general expression 
for the gluon Green's function 



J'u^ifyJ^iito) + S{u})J^^{t)J^^{to) 



t>to. (8) 



G{Y;t,to) = l^e-y'^ 
J 2m uo 

^This form is a standard one for a solution based on the 7-representation (Mellin transform), and 
can be shown to be sufficiently general with the approach of the w-expansion in which K is taken 
to be dependent on the variable w, conjugated to Y . 
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This expression is interpreted as a decomposition into perturbative and nonperturbative 
part of the gluon Green's function 



G{Y;t,to) = Gp,rt{Y;t,to) + G^{Y;t,to) . (9) 

The second term on the right hand side of Eqs. (|]) and defines the Pomeron con- 
tribution at y ^ oo 

Gp(r;t,to) ^ RexpiujjpY) J^^^it) J^^M— ^ (10) 

UJ-p 

where we set R = Res^^jpS^uj). For instance, in the case of the Airy diffusion model 
with cutoff at t = t, i? is roughly estimated [|13[ to be equal R = Aci;p/^o[Ai'(— ^o)]^; 
on the basis of the bound state condition Ai(— ,^0) = 0, with Atup = uj^it) — cjp 
^ [^/DbaM^^'^JP, where D = x'W^Xm- 

Result (^) suggests an estimate of the Pomeron contribution on the basis of the 7 
representation of the regular solution 

where we have used 

a,(t) = i, 6=11/12. (12) 

This representation was argued to be valid in the perturbative regime 
for kernels of type for which ^'(7) = xil) is the characteristic function of K, and 
can be extended to the case of more general kernels in the framework of ci;-expansion. 
Suppose now, that besides ha sit) ~ bu-p also is a small parameter. This is possible in 
a realistic case (that is b ~ 0{1)), if, for some reason tup itself is small (weak Pomeron) 
and is also possible in the case of fixed cup in the limit of 6 — > 0, i.e. in the small 
6-expansion being considered here. Then, one can estimate the magnitude of JF^^ from 
the saddle point condition applied to (0) 

buj-spt = xi^-p) (13) 
which implicitly defines p{bujpt) and yields 



jr^p(t) ~ exp ( / p{x)dx] J . ,M , Xm = x(l/2), (14) 



buJTP / V -27rx'[p(6t^pt)] 

where we have used the relation 



p{bujpt) bt H X{- - p{buj-pt)) = — / p{x)dx . (15) 
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Because of the decomposition and of the perturbative behaviour (|^) we obtain 
the estimate 



btR 



Gpert{Y]t,t) ~ \x'{l)\uJwS TT 



exp 



buj 



p 



p{x)dx 



■ (16) 



In other words, the Pomeron is potentially leading for Y ^ oo because tup > ujs{t) 
for t ^ 1, but is exponentially suppressed with the hard scale parameter t. We notice 
that the suppression exponent is dependent on the parameter buj-pt = cup/a^, and also 
on the model being considered (implying a different solution for p{x) from (p!3|)), but is 
always of the form 



bas{t) 



;i7) 



(the uj-p dependence is implicit in g). For instance in the case of the Airy diffusion 
model we get 



X = Xm[l + Dp{x) 



p{x) 



X 



/D V Xm 

Therefore in the limit bu^t 1 we get the suppression factor 

4 



- 1 



;i8) 



exp 



UJjp 



Sbos Y DuJsif) 



exp [—const 



(19) 



with the function g ~ const/y^. On the other hand in the collinear model with 

xil) = 1/(7(1 ~ 7)) obtains 



1 



X 



i -p(x)2 

which yields the suppression factor of the form 

1 



p{x) 



1 



(20) 



exp[— -— ] ~ exp(— t), g ~ const . 



(21) 



We conclude that, at large values of t, the transition probability to the Pomeron 
behaviour is typical of a tunneling phenomenon and is of the general form 



exp 



ba 



(22) 



so that the b parameter at fixed values of plays a role of /?- in a semiclassical approach 
in quantum mechanics. 

A better insight into the tunneling phenomenon is gained in the simple example of 
the Airy diffusion model which is obtained by expanding the kernel eigenvalue into the 
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second order in 7. By using this approximation, the evolution equation (|1|) takes the 
form 



G{Y;t,to) = iJsit) 



1 + D 



G{Y;t,to). 



We make then the following change of variables 

"* dt 



X 



which in this case results in 



X 



2 



(23) 



(24) 



(25) 



Making use of (|2^j25|) and performing additionally the transformation G = x eG one 
can recast Eq. (ESf) into the Schroedinger-like diffusion equation 



^ G 



+ V{x) 



G 



(26) 



where now the rapidity Y is interpreted as the imaginary time and p 
"momentum" variable conjugated to x. The potential V{x) has the following form 



ibS- is the 

ax 



V(x) = -LJ,[t(x)] 



36 



x^ 



(27) 



and is basically provided by the —LJs[t{x)] part dependent on the running coupling as{t), 
apart from the quantum corrections proportional to b^. 

In this interpretation, which differs slightly from the cj-dependent ones discussed 



previously |]I2| , the non-perturbative parameter is directly given by the lowest energy 
state in the stationary equation corresponding to (p6D, whose potential is pictured in 
Fig. 0. Therefore the barrier is directly that arising at energy — a;p by the fact that 
ujs{t) decreases to zero at large values of t{x). Both the binding of the Pomeron and 
its suppression at large values of t are provided by the running of as{t) which implies 
that the strong coupling region is much more "attractive" than the weak coupling one. 
The detailed estimate of ujp depends crucially on the shape of the potential well for 
small values of t which instead is provided by specifying the strong coupling behaviour 
of ti;s(^) and ds{t). For instance in the case of the cut-off prescription (Fig. P 



asit) = e(t - t) 



bt 



one requires the solution to vanish at the boundary t, 



(28) 



(29) 



The suppression factor on the other hand is basically provided by the perturbative 
regime and by Eq. ( [T^ ) in the semiclassical approximation. 
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Figure 1: Shape of the potential V{x) in Schroedinger-like equation from the diffusion 
model with running coupling, Eq. Solid line: leading term —Usit{x)); dotted line: 

leading term + quantum corrections — The dashed line illustrates the boundary 

x{t) at which one has specified the regularisation of ds{t). 

In the potential-like interpretation above, the Pomeron contribution to the Green's 
function corresponds to the lowest lying bound state and nearby ones (which are so- 
lutions to the homogenous equation) and the perturbative contribution roughly cor- 
responds to the continuum of positive energies, yielding singularities at Rea; < 0, see 
Fig.H. This decomposition is however not unique (neither was (|^), because of ambigu- 
ities in the definition of the irregular solution JF^). In fact, subleading bound states, 
Uj, < Up, see Fig. ^ can be incorporated into the perturbative part also, depending 
on the size of Ugit) . Since they are still suppressed by exp(— ^^i^) factors, like the 
Pomeron, they appear to play a role similar to renormalon ambiguities in QCD pertur- 
bation theory. The important point is that the whole Green's function is well defined, 
given the strong coupling boundary condition, as is also the Pomeron contribution, 
provided by the lowest energy bound state —ujjp in the potential V{x). 



9 



-co 




Figure 2: Schematic representation of the states in an arbitrary potential well. The 
continuum, which corresponds to Re uj < 0, is the perturbative part; is the lowest 
lying bound state corresponding to the Pomeron; Uf, denotes the subleading (discrete or 
continuum) bound states. 

3 Perturbative Green's function and its 6-expansion 



We have just understood on the basis of a few solvable examples, that the gluon Green's 
function can be decomposed (in a non- unique way) into a perturbative part, which gives 
the uJs{t) exponent, Eq. (|^) at large values of t, and a Pomeron part which is leading 
for F — > oo, but is suppressed by the tunneling factor (|1^) at large t's. We shall now 
analyse a similar decomposition in the more general case of a BFKL model with running 
coupling, in which K has the form (|^) where is a scale invariant kernel with eigenvalue 
function xil)- K is known |]10[ that for such a model, a formal solution for the Green's 



function in the perturbative regime is provided by the double 7-representation 



yc,{t,to) = — exp 

_ 2m .L, UJ 



7 



(7 - \)t - (70 - \)t, - ^^(7) + ^^(70) 



(30) 



This expression should be used with care for several reasons. Firstly it is a formal 
solution of Eq. (|l|) with running coupling as{t) = 1/bt which diverges at t = (Landau 
pole). Thus it may be acceptable for t.tQ^t only, where t parametrizes the boundary 
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of the strong couphng region, when as{t) departs from the perturbative expression 1/bt. 
Furthermore, the representation ( pO]) has an essential singularity at = 0. Therefore, 
the Pomeron singularity cjp > arising from the regularisation of the strong coupling is 
not obviously included in this representation. Finally, the definition (^) is ambiguous 
also, because the convergence path of the 70 contour is not specified for 70 — >■ cxd. Any 
further specification will represent the addition of a single 7- representation ([TT|), that is 
a solution of the homogenous equation (|l]) (as is the Pomeron). Given these questions 
we shall first provide a more rigorous definition of the perturbative part of the gluon 
Green's function. 



3.1 Relation of the double-7 representation to the spectrum 

We start considering the double 7 representation ( pUD in the context of the principal 
value regularisation of the Landau pole in Eq. (^. Such regularisation was already 
considered in Ref. [§] and was shown to result in a continuous u spectrum, with 00 
having all real values from — cxd to +00. This seems paradoxical at first sight, but is 
actually hardly surprising, since in such a case the strong coupling constant as(t) is not 
positive definite and is unbounded. Since we actually aim at a physical Landau pole 
regularisation in which as(t) keeps its positivity and is bounded, we should properly 
handle in Eq. ( pO| ) the contribution of the positive end of the tu-spectrum, which is an 
artifact of the principal value regularisation. 

Let us recall ^ that, in the analytical principal value regularisation^] 

(31) 



t J 2 \t + ie t — ie 
the eigenfunctions in 7 = | + zi/ space of the kernel (^ take the form 

Ml) = f^M-^xi',)], 

where we have used the representation t = —d^ to derive the (right) eigenfunctions, 
and the conjugate (left) eigenfunctions due to the asymmetrical role of the coupling in 
Eq. (^, are defined by ^ujilo) = ~^-^i2(7o) ^.nd satisfy the delta- function normalisa- 
tion [T^^T^i) = 6{uj — uj'). As a consequence, the Green's function = {u — K)^^ 
satisfies the spectral representation 

QUl^lo) = / j^u;'{l)^u;'{lo) , (33) 



^ Other prescriptions may lead to different results. For instance taking a sharp cutoff 
p(t) ^ \Q{\t\ - e) separates the Hilbert space into the positive-t and negative-i eigenfunctions, with 

a distortion which, for large t values, is again O ^exp(— — 
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where the variables 7 = 1/2 + zz/ and 70 = 1/2 + iuQ run along the Re 7 = 1/2 axis. In 
order to evaluate the integral (|33|) we shall distort the u' contour so as to go to Im cj' > 
(Im u' < 0) for z/q > (z/ > z/q). We have to additionally assume that Xm < 0, so that 
z(X(7o) — X(7)) has always the sign of z/q — u. The case Xm > will then be treated by 
analytical continuation. The result of this contour integration in Eq. (|33D reads then 



^^(7,70) 



1 d 



ec^(t'o - t^) exp 



(X(7o)-X(7)) 



with 




Im uj > 
Im cij < 



(34) 



(35) 



Then, by performing a Fourier transform to t-space, it is straightforward to derive 
from the expression (|3^ ) the double 7-representation of Eq. ( |38D with the additional 
constraint that the contour c(oo) is chosen so as to make Re (70 — 7)/^; > 0. 

It is now important to realise that the positive u spectrum in (p3|) is strongly 



dependent on the details of the regularisation procedure around the Landau pole of 
the running coupling as{t). In fact, if one assumes that Xm < and as(t) is frozen 
for t < t where i is large and negative, one realises that the spectrum has a gap for 
< LJ < Xm/ipt). The simplest way to see it is by using the Schroedinger picture of 
Ref. [^, in which the eigenf unctions are solutions with negative energy Xm in the linear 
potential]^ ~ hojt plotted in Fig.^. We need in fact that the potential depth \bujt\ > \xm\ 
for the continuum to begin. Because of the gap, the Green function ^a;(t, to) can 
be continued analytically in the whole cut tu-plane, while the eigenfunctions and the 
expression ( ^^ will be only slightly distorted for t, to > if — t 3> 1. Therefore we can 
define, see Fig.H 



GpertiX'-, t, to) 



•Y 



(36) 



This expression cuts off the (Pomeron-like) positive u spectrum completely and 
yields an asymptotically decreasing function for both positive and negative Y^s (recall 
that Xm < 0). Furthermore, for |t| ^ 1, one can use (|34[) and therefore the double 
7 representation (^) in order to estimate it. We have thus achieved our goal, that 
was to show that the double 7 representation, with the tu-contour running over the 
imaginary axis, is an acceptable definition of the perturbative part Gpert{Y;t,tQ). In 
fact, while the negative u spectrum with large eigenfunctions in the large t ^ 1 regime 
is weakly dependent on the regularisation procedure, the positive cu spectrum strongly 
depends on it but it is completely cut off by the uj integral being used. Since the a; > 
eigenfunctions are large in the t < regime, changing the regularisation procedure to 

■^Strictly speaking, this picture applies to eigenfunctions in the Airy regime of large t, with 
{bujY^'^{t — ^) kept fixed. That's enough to hint the spectral properties we need. 
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Vo,(t) 




Figure 3: Plot of the linear potential boot with freezing below i. Xm has here the inter- 
pretation of the energy in the Schroedinger-like problem. Dashed lines correspond to the 
principle value regularisation of the running coupling. 

a physical one will now affect the perturbative part only through the tunneling from 
t < to t > . Such effect is suppressed as ~ exp[— l/(6Q;s(t))] for \t\ ^ 1, as was the 
Pomeron in the estimate in Sec. 2. 

In other words, we propose to define the perturbative Green's function by the u- 
contour over the imaginary axis in Eq. ( PB] ) in the double 7 representation with Xm < 
(followed by the analytic continuation to Xm > 0). This is an acceptable definition and 
cuts off the positive o^-spectrum in the principal value regularisation, to be replaced 
by the Pomeron one in the physical case. However there is an intrinsic, renormalon- 
like ambiguity in such a definition due to the fact that the strong coupling boundary 
condition still distorts the perturbative eigenfunctions by terms which are of the relative 
order ~ exp[—l/{bas(t))]. 

3.2 6-expansion and diffusion corrections 

We have just shown that the gluon Green's function for the kernel with a physical 
regularisation of ctsit), can be written as, 

G{Y;t,to) = Gpert{Y;t,to) + Gp(F;t,to), (37) 
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CO 



Figure 4: Spectrum of uo values with a gap and the corresponding contour for the per- 
turbative Green's function. 

where we propose to take 



(38) 

and c(oo) is a proper specification of the 5 contour with Re- > . As was mentioned in 
the previous subsection this decomposition is not unique, but here we focus on properties 
of the perturbative expression (^) which are independent of renormalon ambiguities. 

Let us now evaluate ( |38D in the regime in which t, to, ojY are all large param- 
eters. One way of looking at it is to let 6 ^ with as{t),as{tQ) and hY kept fixed, 
which is the 6-expansion we are investigating. This allows us to perform a saddle point 
estimate of the exponent in (|38|) for any uj value. It also ensures that the Pomeron 
contribution is strongly suppressed by the tunneling factoiQ (0)- In such a regime we 
obtain the saddle point conditions 

hut = xil) 

biJto = xil + S) (39) 
biJ^Y = X(7 + 6)- X(7) . 

For instance, if we take to = t and we use the (anti)symmetry of xil) (^(7)) fo^' 
7 ^ 1 — 7, we can define l = \ — p with the equations 



5 = 2p 

2 



buj'Y = -2X(i-p) = 2X(^+p), (40) 



'^In particular regimes of the Y and t parameters the present perturbative estimates are still valid 
for b — 0(1) but the Pomeron suppression is not strong and contaminates the perturbative behaviour. 
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Figure 5: Plot of the g{p) function ^j^ ) in the case of collinear model xil) 



7(1-7) ■ 



and the saddle point exponent of (|3^) becomes 

2 1 2 I'^^i 

E{Y;t,to) =LuY -2pt + —X{-+p) = ujY - — p{x)dx , (41) 

to 2 bu J^^ 

where we have again used the integration by parts, analogous to (|T5|). Insight into the 
behaviour of (^Tj) is obtained by eliminating lj in Eq. to get (uj'^it) = Xm/bt ) 



(jj = 



bt ' 

bt^ xni-p) 

We see, Fig. |], that for any given value of parameter ( = u'^{t)Y/t = bagUj'^Y fixed by 
both energy and scale, there is a value of p which in turn determines the saddle point 
uj = as{t)x{^ ~P[C])- When ^ <^ 1 is a small parameter, the solution of (^21) is driven 
towards the minimum of x ? corresponding to p = 0. The other solution with p ~ 1/2 
(7 ^ 0) is instead subleading and unstable and is thus discarded. By expanding the 
eigenvalue function x around the minimum we obtain, 

1 _L 1 n^2 p. „ „ (4) 

^ ^ I n 2^2 ^ 2p(l - '-Df + {by + 
(1 + Dp^y 3 3 5! Xm 

p ^ ^(1 + 1/^^2 3 1 

^ 2^ 12 ^ 8 ^ SO 4\xm 

/•2 c 1 (4) 

. ^ .«W(l + 44flV + i|^C«...-), (43) 
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and the saddle point exponent (^T|) becomes 



E{Y;t,t) 



1 1 1 v^^^ 

12 ^ 24 ^ 120 Xm 



(44) 



The result exhibits the existence of a power series in the parameter z = D^"^, 
previously emphasised JTSf in the context of the Airy model. The first correction of 
this series yields the term ~ D{uj'^Y)^/t'^ found in various papers, for example see 
0, H, 0. Further corrections, subleading in Y, come from integrating the fluctuations 
around the saddle point (|40|) . The fluctuation matrix has two positive and one negative 
eigenvalue, roughly corresponding to 7 and uj fluctuations along the imaginary axis and 
to 6 fluctuations along the real axis (see Appendix A). As a result one obtains the 
following normalisation factor 



Nit,Y) 



-47rx'(7) 



1 + 



X'(7), 



Xm 



'l + aDe + ---) 



with 



5 

12 



1 (4) 

24 (x:;)^ 



(45) 



which shows further dependence on z = -D^^, of the type ~ Y^. From higher order 
fluctuations one finds also further subleading corrections with linear dependence on Y. 
In appendix B we show how to obtain them in the specific case of the Airy diffusion 
model (see Sec. 4). In the BFKL case the term linear in Y is calculated in and 
provides the shift 



12 



15 

T 



XmX 

4 



(4) 



"0 

m 



1 XmXm ^ I XmX^ 
4 -^"4 8 A/"3 



Xr. 



Xr. 



(46) 



0{1) and < 
and to the collinear behaviour 



to the hard Pomeron exponent due to the diffusion corrections. 

Let us note that the saddle point (^21) can also be studied for ( 
p < |. Due to the expansion of x around 7 = 
^(7) — ^ + + ■ ■ ■ ) the function g{p) in the r.h.s of (|42|) is nonnegative and vanishes 
for both p — i> {g{p) ~ cp) and for p — > |, ( g{p) ~ (p — l/2)^log(p — 1/2)), see 
Fig. ^ Therefore it has a maximum at some value of p = pc- Correspondingly, for 
( > (c = gipcj , the saddle point(s) become complex conjugate, and the exponent 
E{ujY, () acquires a branch point. In the case of BFKL the critical value Cc — 0.264. 



Such behaviour was hinted at in ||T8| and is here shown to be quite general. 

To summarize, we have shown that the perturbative Green's function in the b- 
expansion takes the form 

n{C,t/to) 



GpertiX'i t, to) 



2^/MJf{t)Y 



exp[a;°(t)r£:(C,tAo)] 



(47) 
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Figure 6: Upper plot: the saddle point solution u{t) as a function oft for the collinear 
model. Lower plot: saddle point solutions jlt) (dashed line) and '^^{t) (solid line) . The 
parameters are : b = 1, rapidity Y = 5 and scale to = 15. 

where C = Ug{to)Y/to and for t = to, 

n{0 = l + aDC + ..., s^Q = i + 1dC + ... , (48) 

are given by power series in ( with convergence radius Cc which is equal to the maximum 
of the function g{p), see Fig. |^. Note that the exponent in Eq. (^Tf ) is O (^) at fixed 
values of bV, as{t) and (, while the normalisation is O {1) and the fluctuations are 
O {b'^), n > 1, as is natural from the role of 6 as a semiclassical expansion parameter. 

The above presented analysis can be easily extended to the case of different scales 
t > to, in order to reproduce the collinear limit. The exact solution to the set of 
equations ( |39D is rather complicated to be found analytically for the general case of 
non-equal scales; it can be however easily solved numerically. In Fig. |^ we show the 
solution for the lj, 7 and 70 = 7 + 5 in the case of the collinear model xil) = V(7(l~7)) 
as a function of scale t starting from t = to- The saddle point exponent 00 first decreases 
with t, then has a minimum ujm = ujf'\to) at certain tm and then increases. The solution 
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Figure 7: Function (7(70; t, to) from Eq. ( f5j| j plotted versus 70 for different values oft: 
15,15.2,15.5,16,18,20. The value of the maximum increases with increasing t. The 
dotted line represents the value of the parameter Q = XmY/ (btl) . The parameters are 
the same as in Fig. |^. 

for the 7o starts from the value 70 > 1/2 but then moves towards zero as the 7 also does. 
The point at which 7o(tm) = 1/2 is exactly the minimum of the saddle point exponent 
uj(t). In the regime when t to ^ 1; both solutions tend to zero (7,70 0) and one 
can recover the correct asymptotic behaviour by making the coUinear approximation in 

m 

1,-1 

bout ~ — , buto ~ — 
7 7o 

bu^Y ~ln^~ln- (49) 
1 to ^ ^ 

which gives 



a;^^/^ln|- and E(y; t, to) ^ 2^1^ In |- . (50) 
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The function g{p) is generalised in this case too, and depends additionally on t and 

to 



9{lo\ t, to) = Xm7o(l - 7o)^ In 



7o 



I + a/^a/(7o - 1/2)^ + f. 



1 - 7o 1 

2 



(70 - 1/2)2 + A,) 



(51) 



where we have defined At = t — to- In Fig- we show it as a function of 70 for different 
values of the scales t and to- Again the value of the parameter ( = Xmy/{btl) = 
5f(7o, t, to) will determine the solutions for the 70 and the related quantities. The solution 
will be given by the lower value of 70 and it will quickly pass 1/2 and shift towards 
with increasing asymmetry At. The maximum of the curves (7(70; t, to) also increases 
with At which means that in the coUinear situation when t ^ to ^ 1 one has a real 
solution for quite large positive ( and that the critical value Cc ~ Int/to. Note, however 
that due to the pecuhar shape of the function g{'jo;t,to) (Fig. 0) for ^ — 0"*" both 7 
and 7o drift towards zero, as in Eq. (ff9|), and that the exponent in Eq. (pO) implies a 



singular behaviour of the exponent function S in Eq. (47) 



S{C,t/to) 



'logt/to 

XmC 



C = o;,(t)F/to 



(52) 



Therefore, for t > to, the exponent function S is no longer a power series in ( around 
C = 0, but should be rather expanded around some nontrivial value of ( (like ( = 
g{^; t, to), for which 70 ~ i), so as to stay away from both C = and ( = Cc- 

On the other hand, if (t — to)/to = bas{to)At -C 1 the exponent in Eq. (|50|) takes 
up the double log DGLAP form ~ 2^yYas{to)At with frozen coupling (one cannot 
distinguish as{t) from as{to) for such values of At). Since this expression is finite at 
fixed as{to)Y, one can still perform an expansion of the exponent in b in this small At 
regime, as proposed for a numerical procedure in Sec. 5. 



4 Explicit Green's function in the diffusion model 

The diffusion model with running coupling of Eq. (|23| ) has been discussed in the liter- 
ature |2^ in various ways, starting from its equivalence |29| with a Schroedinger-like 
problem in t space. We reconsider it here with a purpose of providing a representation 
for it which makes the 6-expansion and the evaluation of diffusion corrections quite 
explicit. 

We start from Eq. (|38|) , with Xm < and with the specification of the 6 contour 
provided by Eq. (|3^), according to which Re{6/uj) should be positive. The exponent in 
([38| ) in the case of the diffusion model has the following form 

5y 1 8-y" 1 1 y" 

E:=ujY + {^- l/2)At - 5to + ^ + -^(7 - 1/2)^ + -^(7 " 1/2) + -^5' - 

bu 2, buj 2 bu o buj 

(53) 



19 



We notice that, since XIj) is cubic in 7 for the diffusion model, X{'-f + 6) — X{'-f) 
is quadratic in 7, with coefficient of the quadratic term given by which has a 

positive real part. Therefore the 7 integral, at fixed u and S/uj, converges along the 
imaginary axis and can be done explicitly to yield 



Gpert{y\ t, to) = / — — / dx 



btr 



m 



exp 



2/^ + ^0 X , 

- X (u— aoXm) + 



2t 







(54) 



where we have introduced the variables 



X 



At = t-tn 



UJ 



1 y" 

Am 



6 bto4 



2 ' 



"0 



1 

bk 



(55) 



Surprisingly enough, the uj dependence of the exponent in Eq. (Q) is still quadratic at 
fixed ~ 5/ti;, with positive coefficient of the quadratic term. Therefore the uj integral 
converges along the imaginary axis and can be done explicitly to yield 



G pertly '1 



dx 



btr 



exp 



(Y-x^^-±tr 



n{to)x^ 



~bt^ 



{Atfbto 
2v" 



(56) 



where we recall that the x integral converges because we have set Xm < . Finally, by 
introducing the new integration variable ^ = we obtain (t = ^^) 



^0 f' 
Gpert{y',t,to) = y / 
t Jo 



X exp 



v/2vrx'4a,(t)|F| v/vrfi(t)|y| 



{Atn 



n{t)\Y\ 



^X'Lo^s{t)\Y\ 



(57) 



where Ugit) = ^ < and the — (+) sign holds according to whether F > (F < 0). 



Several features of Eq. (|57|) are worth noting. First, the solution decreases for 
\Y\ — ^ cxo in both directions, as expected from Eq. (|5BD, but more strongly for Y —00. 
Secondly, it satisfies the following boundary condition 



G{Y = 0+; t, to) - G{Y = 0' ; t, to) = S{t - to) , 
but the Y = 0~ contribution is non-vanishing. As a matter of fact we have 

POO 

G(F = 0+;t,to) = / rfe5(v^(t-to))5(v^(e-l)) = 25(t-to) 

POO 

G(F = 0-;t,to) = / d^5i^it-to))6i^i^ + l)) = 6it~to) , 

Jo 



(58) 



(59) 
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and this means, according to Eq. (p6D, that the projector over the negative (positive) 
spectrum is 2(-l) in this case. Third, apart from the overall 1^/1 factor(due to the 
asymmetrical role of the running coupling in the kernel), the t and to dependence occurs 
mostly through the parameter t = corresponding to the scale fcfco for the hard 
process. Finally, the 6-dependence is nicely summarized by introducing the parameter 



(60) 



in terms of which Eq. 



GpertiX] t, to) 



reads 



i Jo 2ny/D\z\/3 



exp \ uJs(t)\Y\ 



i-z)/3 



Aty ^ 1 



(61) 



for Y > {Y < 0) respectively. This equation provides the almost explicit representa- 
tion of G we were looking for. 

Before analyzing Eq. (pl| ) in more detail, let us consider the analytic continuation 
of Eq. (|57D to the positive values of Xm- Only the ^ = region of the integral is 
affected. In this case the contour in the ^-complex plane has to be distorted so as to 
reach the end point ^ = from the Re^ < region, and this can be done in several 
ways, see Fig.|^. Since the measure of the small ^ region vanishes, this contribution is 
not expected to introduce sizeable ambiguities, which are roughly given by solutions of 
the homogeneous equation with u = 0. 

The large positive Y behaviour of Eq. (pl| ) at fixed value of z is determined by the 
saddle point 



/(e 



At, 



(3^ - 1)3(^ - 1) 



At 
2^ 



2n 



(62) 



For z < (xm < 0) there are four saddle points, the leading one being close to ^ = 1. 
For z > (xm > 0) the saddle points close to ^ = (contributing to the completeness 
relation) become complex and two real ones remain. The leading one at ~ 1 — e for 
At ^ t and 2; -C 1 is equal to 



1 



4i) + + 



At 
2^ 



(63) 



and yields the exponent 

E{Y;t,to) = uJs{i)Y 



1 + 



z{i) , z{if 



12 



+ 



24 



At 
"2! 



zit) 



(64) 



thus confirming the results of Ref. [|T8| and generalizing them to arbitrary scale depen- 
dence. 

By taking into account the fluctuations around the saddle point we can also obtain 
the ~ Y^ corrections (apart from the already known ~ Y^ ones of Eq. (p^). It turns 
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Figure 8: A possible continuation to Xm > of the contour in the ^ plane for the integral 



( pip representing the perturbative solution of the Green's function in the diffusion model. 



out that by considering the third order derivative of the exponent in (|61|) one can also 
identify the correction hnear in ~ F which corresponds to the second order shift of 
ujs{t). The details of the calculation are provided in appendix B and yield the result 
asXnX, which reproduces the first term of 

The importance of the parameter z{t) is confirmed by the above procedure. In fact, 
above some critical value Zc{At/t) the saddle points collide and become complex, imply- 
ing a singularity of the saddle point evaluation and thus a change of regime occurring 
for Us{t)Y ~ V^- '^^^ critical value of Zc is ^-dependent and ranges from |(| + 
for At = to = for ^ = 1, (to <^ t). This has to be contrasted with the results 
for the case of collinear model (see subsection 3.2 and discussion about Eq. (plD) where 
the critical value of z^ is always positive and has the asymptotic behaviour z^. ~ In^t/to 
for t ^ to ^ 1- The fact that z^ = for the Airy model in the collinear regime just 
means that this approximation is no longer valid. In fact the validity of Airy model is 
strongly limited to small values of At . Finally, for z ^ 1, the complex saddle points 
at ^ = ±^(e^*'^ z)^/^ dominate the exponents, which become in that case 

EAY;t,to) = -i...(t)y.-^e^x = -^y^^e-^-f , (65) 

thus implying an unphysical oscillatory behaviour for uJs{t)Y ^ t, as already discussed 
in m. 
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Figure 9: Plot of the "saddle point" function f{C,, At/t = 0) versus C,, The dotted 

line represents the value of —z variable, ^^). The intersection points define the saddle 
point solutions for the exponent in solution (1^7[j. 

5 Numerical results 

In addition to the analytical studies carried out above, it is also possible to examine 
the 6-expansion from a numerical point of view. There are two main purposes to this. 
One is to examine the structure of the 6-expansion somewhat more generally, be it at 
higher order, or in the context of more general kernels. 

The other aim is to establish a way of numerically defining 'purely perturbative' 
predictions for high-energy scattering, as well as the potential domain of validity of these 
predictions. While for leading-order BFKL calculations this is not strictly necessary 
given the analytical tools at our disposition, when including higher-order corrections, 
numerical methods may represent a more practical approach. 
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5.1 Numerical 6-expansion 



The BFKL equation (|l]) for the gluon Green's function can be written in rapidity and 
transverse momentum space as 



dG{Y;t,to) f d'q 



dY J 7iq 



-.2 



^G{Y-Mk'\t,)-G{Y-tMmk~q) 



k' = \k + q\, 
(66) 



where the scale X of the couphng may be t (as has been the case so far in this article) or 
Ing^, or some other combination of scales in the problem. The gluon Green's function 
is defined with the following initial condition: 

G{0;t,to) = 6{t-to). (67) 

To determine the 6-expansion of G, we carry out several numerical evolutions of the 
initial condition, each time with a different (small) value of b. For arbitrary b, the 
coupling is defined such that agito) is independent of b: 

l + (t-to)to.(to) ^ 

Typically we use b values of —n6b, — {n — l)6b, . . . , n6b with 6b of order 0.01 and n = 3. 
The smallness of 5b ensures the absence of non-perturbative tunneling contributions, 
since they are suppressed by a factor exp[—ji^^^-^g{as)]- The use of negative 6's may 
seem surprising, but actually Eq. (|68D is valid for any sign of b, and the 7-representation 
also, provided one replaces to by a~^{tQ)/b. One then considers the 6-expansion for 
lnG(F;t,to), 

In G{Y; t, to) = ^'^^(^' t - as{to)) • (69) 

where we assume that t — to <^ to- Expression (p9| ) corresponds to the series in 
Eq. (^81) for the normalisation and exponent functions in terms of the parameter 
( = bas{to)i^s{to)Y ■ We note that this numerical 6-expansion differs slightly from the 
analytical 6-expansion discussed earlier. In particular while previously the expansion 
was performed with as(to), as(t) and bY fixed, here it is as(to), t — to and as(to)Y that 
are kept fixed. 

In the approximation that for the small 6- values being considered. In G is well rep- 
resented by a truncation of the series (l69| ) at term i = 2n, it is straightforward to 



determine the ii. Thus with n = 3 we are able to determine terms in the expansion up 
to ig. Formally the error on the determination of the ii from neglected terms is of or- 
der {n6b)'^'"'~^^~'^i2n+i, though it can be larger if there is significant cancellation between 
terms.R 



^In practice numerical rounding errors are also important and contribute at the level ofe/6V, where 
e is the relative machine accuracy. The trade-off between these two sources of error determines the 
optimal value of Sb. 
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In a number of the figures that follow we will actually consider the effective exponent 



(jJes{Y; to, to) = ^777 In G{Y; to, to) 
aV 



(70) 



rather than the Green's function itself — the use of Ues facilitates the identification of 
the F-dependence of certain of the £i, and Ues can itself also be expanded in powers of 
6, 



i=0 



(71) 



Figure |1^ shows the uies,i for evolution with scale t and as{to) = 0.1. This low value of 
Us has been chosen so as to have a value of ujsito) — 0.277, close to that expected for the 
more physical asito) — 0.2 once one takes into account NLL corrections. Odd powers 
of b have zero coefficient, due to the b ^ —b symmetry of Eq. (0), which in turn is 
due to the 7^-^1 — 7 symmetry of the BFKL eigenvalue function and eigenfunctions, 
Eq. m. 
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Figure 10: Left: terms in the 6-expansion of UesiX'jtoyto) up to order b^, and addition- 
ally the asymptotic analytical prediction for the b^ term {b'^a^Y'^ and b'^a'^Y contribu- 
tions). Right: the difference between the full numerical result for the b^ term and the 
asymptotic analytical result at lower and higher numerical precisions. In both cases 
asito) = 0.1. 



The O (6^) term in the left-hand plot illustrates the characteristic Y'^ dependence 
expected from Eq. (|44|) (due to the derivative the contribution b'^a^Y^ to InG gives a 
contribution proportional to Y'^ in a;efr,2)- It is compared to the sum of the analytically 
determined b'^a^Y'^ and b'^a^Y terms and asymptotically one sees a good agreement. 
The right-hand plot shows the difference between the full numerical and partial ana- 
lytical determination of ujes,2 and one sees that it is consistent asymptotically with a 
constant term, which equals to 5uj^p /b"^ ~ —0.002491 from Eq. (^). The two curves in 
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the right-hand plot have been determined with numerical evolutions of different accu- 
racies so as to illustrate the numerical stability of the procedure. The differences visible 
at smaller Y values are a consequence of sensitivity to the different approximations 
(widths) of the initial delta-function in t. 

The left-hand plot also shows the ci;eff,4 and ujesfi terms. While at low Y they are 
suppressed, they grow much faster with Y (as Y^ and Y^ respectively) and quickly 
dominate over the ujes,2 term. Of course, they are suppressed in Eq. ( |7TD by a power 
of b also. But if we take b = O {1), then the ( parameter in Eq. (|4^) is sizeable, and 
the importance of higher order terms increases with Y, meaning that the series (|7l|) is 
slowly converging, because of nonanalyticity at C = Cc- This point will be discussed in 
more detail below and determines the perturbatively accessible range of Y for BFKL 
predictions. 
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Figure 11: The difference between uJes^i{Y; Iq, to) for the case of a coupling that runs as 
^^(Ing^) and one that runs as as{t). Shown for as{to) = 0.1. 



It is also of interest to examine what happens to the evolution when as is evaluated 
not at scale t but at Ing^, which is favoured by the NLL corrections to the BFKL equa- 
tion 1^. The differences in the Ues,i between these two options are shown in figure |TT]. 
This figure has been generated with a Gaussian (oc e~*-*~*°^ ^'^) initial condition as op- 
posed to a delta-function, because with a delta-function, at small Y one finds spuriously 
large contributions enhanced by logarithms of the 'width' of the delta-function. 

With Q;s(lng^) one loses the b ^ —b symmetry of the evolution and both odd and 
even powers of b are present. At orders b and 6^ the effect of changing from as{t) to 
as(lng^) is simply to modify Ues by a constant — essentially the dynamics which led 
to the asymptotic F-dependence in figure 1^ is independent of what scale one chooses 
for tts- The constant shift is trivially A = —u^^ at order b, and contains both scale 
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changing and diffusion effects at order 6^. It is only at order and beyond that Y- 
dependence starts to appear, and one sees a form of mixing between the shift in Ues 
which appears at relative order bu^ due to the scale ^^(Ing^) and the dynamical effects 
which arise for any choice of scale at relative order (feagF)^. This is reflected in the fact 
that the order b^ and 6^ terms have a leading dependence, while the b^ and b^ terms 
have a leading ^ dependence. 



5.2 Pure perturbative predictions 



One application of the 6-expansion is that of extracting 'purely perturbative' numerical 
predictions. 

One of the most common ways of obtaining BFKL predictions including running 
coupling is to solve Eq. ( |66D (or its analogues with various forms of higher-order correc- 
tion) with a regularised running coupling, see for example |^6|, 27| . A measure of the 



perturbative uncertainty on the prediction can then be made by comparing different 
regularisations and seeing how they affect the Green's function. 



This is illustrated in Fig. |T2|a, which shows the Green's function evaluated by direct 
numerical solution of Eq. (^) with two different regularisations of the coupling: in one 
the coupling is set to zero for t < t where as{t) = 0.5, while in the other t is defined 
by C(s{i) = 0.25. In this and other plots in this section, we use as(lng^) and the initial 
condition is a Gaussian rather than a delta-function so as to have sensible behaviour 
for G at small Y. 
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Figure 12: (a) the BFKL Green's function evaluated with two regularisation of «<(, a 
cutoff when ols = 0.50 and a cutoff when = 0.25; (b) the Green's function evaluated 
by truncating the sum Eq. (pQl) at various orders. 



At small and moderate Y there is good agreement between the two curves, but 
beyond a certain point they no longer coincide, and one may deem this to be the limit 
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of the perturbative prediction. However we could quite easily have chosen a different 
set of regularisations to compare to (e.g. with a coupling that freezes below t) and 
one would have come to a different conclusion about the point where non-perturbative 
effects become important; it is difficult to make a strong case for one regularisation 
scheme as opposed to another. 

A second point is that for small values of as{to), in the numerical calculation the 
transition to the non-perturbative regime comes about by tunneling |]T^, |T^, which 
places a limit on the maximum accessible rapidity, Vmax, that scales as l/as{to)- How- 
ever there are arguments that suggest that in the real world (as opposed to a numerical 
solution of a linear BFKL equation) effects such as unitarity |]30| , |3T[| , or the fact that 
the Pomeron is soft, will eliminate or suppress tunneling. In such a situation the true 
non-perturbative limit on the maximum rapidity is believed to scale as l/a1{to) [18|. 
There is however a problem of how practically to calculate a Green's function beyond 
the tunneling point. 

The 6-expansion offers a solution to both these problems. The question of how to 
determine the rapidity where perturbation theory breaks down reduces to that of estab- 
lishing when the series expansion in b stops converging. This is illustrated in Fig. 
One eliminates in this way the need for a subjective decision about a regularisation of 
the coupling. It is still necessary to decide where to truncate the series, but one has an 
analytical understanding of the properties of the series and the mathematical tools in 
such a case are well-established. 

The problem of tunneling is also eliminated, because a truncated perturbative expan- 
sion in powers of b will not reproduce a non-perturbative e^^^"^ factor.^ The maximum 
perturbative rapidity in the 6-expansion is therefore expected to behave precisely as 
discussed in [0, namely to scale as l/a^(to)- 



These ideas can be tested by considering a series of evolutions for different as(to) 
values. In each case one establishes a maximum accessible rapidity, l^max, in two ways: 

1. by examining when two non-perturbative regularisations lead to answers which 
differ by more than a certain threshold; 

2. by examining when the difference between different truncations of the 6-expansion 
differ by the same threshold. 



In the first case we take the two regularisations used for Fig. |T2|a, while in the second 
case we take truncations at order b^ and b"^. We define the threshold as being when 
I \nGa/Gb\ = 0.2, where Ga,b are the Green's function for the different regularisations 
or truncations. 

The results are shown in Fig. 13. There is a clear linear dependence on aj^{to) 
for the case of the non-perturbative regularisation, a clear sign of tunneling being the 
relevant mechanism. In the 6-expansion Fmax rises much more rapidly, in a manner quite 



®One may wonder whether tunnehng might manifest itself in the series expansion through some form 
of renormalon behaviour — it is difficult to answer this question properly without going to inaccessibly 
high orders in b. As is shown below however, in practice this issue does not seem to arise. 
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Figure 13: The maximum perturbatively accessible value of Y, as a function of 1/0:5(^0), 
determined by comparing different non-perturbative regularisations of as, or different 
truncations of the 6-expansion. 

compatible with a proportionality to a~'^{to). We note that the value of coefficient in 
front of a~'^{to), of the order of 0.1, is in perfect agreement with the prediction based 
on the critical value Cc = 0.26413 found in Sec. 3.2, which yields Y = j^o;~'^{to) — 

Some comments are in order. For the case of nearly leading-log BFKL evolution 
that has been studied here, one could equally well have performed a normal expansion 
in powers of ag for fixed asY apart from the need of the Pomeron suppression. Indeed, 
it turns out that (for t = to) the b^ii are functions only of agY and bas, and one can 
therefore rewrite the expansion for In G as 

lnG = J2ibc^sy\{asY) (72) 

1=0 

with Xi{asY) = li/a\ = {agYY^^Ci + {asYydi + . . . , (i > 1). From the point of view 
of Eq. ( ^71) the highest Y powers in Aj ~ {usYY^^ correspond to the exponent function, 
the next-highest ones (~ (a^F)*) to the normalisation, and lower ones to higher order 
fluctuations. Therefore, the 6-expansion is needed to suppress the Pomeron, but is 
actually in one-to-one correspondence with an expansion in at fixed asY . 

For next-to-leading-logarithmic (and N"LL) evolution the situation is different. As is 
well known, the series convergence in as is extremely poor, because of several problems 
stemming essentially from large collinearly enhanced terms. To carry out an expansion 
in Us for fixed UsY is therefore almost doomed to failure. The 6-expansion will on the 
other hand be much more stable because at each order in b one will effectively be able to 
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resum coUinearly-enhanced terms a'^{asY)"^, in analogy with what is done in [pO|, p2 



The detailed behaviour is currently under investigation ||19|| . 

It is important keep in mind that after accounting for higher-order effects, the 
numerical values in figures such as Fig. |l^ will be significantly altered, because on the 
one hand the diffusion coefficient is known to be very different at higher orders [Q, 
and the tunneling behaviour is also expected to be numerically substantially modified. 
Nevertheless, we expect the qualitative features of that figure to persist. 

Finally, one potential practical limitation of the b expansion is the case when t and 
to are different, leading to the need for a coUinear resummation. In such a case, as 
mentioned in Sec. 3. 2, in the regime where t — to ^ ^o^ the coefficients of the b expansion 
will be enhanced by powers of {t — to), and though the series may still be formally 
convergent, it is not clear whether it will be practical to include a sufficient number 
of terms for the convergence to be reached. To establish this point requires further 
investigation. 



6 Conclusions 

In this paper we have studied the properties of the gluon Green's function in the case 
of the small- a; evolution equation with running coupling. In general the solution can 
be decomposed into perturbative — hard Pomeron and non-perturbative — Pomeron 
components. The hard Pomeron is then governed by the perturbatively calculable 
saddle point exponent ujs{t) which is modified by the corrections due to running coupling 
effects. The non-perturbative part has instead a true singularity which is leading at 
large values of rapidity since cup > Usit). The Pomeron part is however suppressed at 
large values of t by a tunneling factor which has a universal form of exp[—l/{bas)g{<ys)]- 
As a simple illustration, in the Airy diffusion model the Pomeron singularity is the lowest 
lying energy state of the potential in the Schroedinger-like problem. The perturbative 
part corresponds then to the continuum of states at Re a; < 0. The decomposition into 
perturbative and non-perturbative parts is however non unique in the sense that some 
subleading bound states Uh < ujp can be included in the perturbative part of the gluon 
Green's function and are similar to a renormalon ambiguity. 

Starting from this qualitative picture, we introduced the 6 — >• 0, fixed limit in 
which the Pomeron is suppressed, in order to study the properties of the perturbative 
part Gpert{y',t,to). By applying the b expansion we obtained the saddle point expo- 
nent and the diffusion corrections. We showed that apart from the well known ~ a^Y^ 
corrections there are further subleading ones of type ~ Y"^ and ~ Y. The latter con- 
stitutes a second order shift, which adds up to dynamical corrections of the same 
order b'^a'^Y coming from subleading corrections to the kernel. We have investigated in 
particular a scale change in the running couphng [jl9|. We found out that this pertur- 
bative expansion is controlled by the parameter z = ^{bdsY {ujsYyXm is valid for 
z 1. Outside this regime, when z > 1, the perturbative solution starts to oscillate 
and the non-perturbative Pomeron takes over. Thus the genuine hard Pomeron part 
of the gluon Green's function which is governed by ujs{t) can be studied phenomeno- 
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logically at large values of t and moderate rapidities, in the limited regime where the 
non-perturbative part is strongly suppressed. 

Additional NLL corrections are going to change the picture of diffusion and tun- 
neling in a quantitative way. For example the so called kinematical constraint [p3| , 
incorporated in some resummations of subleading corrections |]T^ can be shown to in- 
crease the rapidity where tunneling occurs by about t — to; the proportionately larger 
effect of higher-order corrections at larger values of will also delay the onset of tun- 
neling, by reducing the non-perturbative Pomeron exponent. Both these effects can be 
expected to widen the window for phenomenological tests of the purely perturbative 
hard Pomeron. It should be also remarked that unitarity effects can in principle change 
significantly the phenomenon of tunneling. It was noticed in [^, |3l| that in the case 
of the non-linear small x evolution equation |3^, 35], the generation of the saturation 
scale Qs{x) leads to the suppression of diffusion into the low scales k < Qs{x) and the 
distribution of the momenta is driven towards the perturbative regime. More detailed 
theoretical studies of these phenomena are thus needed. 
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Appendix A: Calculation of the fluctuations 



The fluctuation matrix is obtained by taking the second derivatives with respect to 7, (5 
and u of the exponent 

E(7, 6, uj)=ujY + {^- l/2)t - (7 + 5 - l/2)to - ^^^(7) + 7^^(7 + 5) (A.l) 

oui boj 

in the double-7 representation Eq. (|38|) . It reads as follows: 

/-2x'(7) -X'il) \ 

-X'il) -X'il) -bt (A.2) 



1 

buj 



\ -bt 2bY j 

where we have made use of the relations (^Ol). The secular equation reads then 
'2x'(7) 



buj 



+ A 



buj J \ ^ 



+ 



buj J \uJ 



x'il) 

buj 



. 2r\ 

A-— = 0. 



(A.3) 
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By looking at the separate terms in Eq. ( |A.3| ) it turns out that the sum and product 
of the eigenvalues have the following form {x'il) < 0) 



Ai + A2 + A3 
A1A2A3 



3x'(7) 2F 
— ^AiZ + _ > 

^[(6t)^ + 6rx'(7)]<0 



(A.4) 



which means that we have two positive eigenvalues and one negative - provided that 

{htf>-^{^)hY . (A.5) 

Since we have hY = ({btY/xm, the condition ( |A.5| ) is equivalent to Xm > —(x'i'l) and 
holds as long as C is a small parameter. Two positive eigenvalues correspond thus to 
fluctuations along the two imaginary axes and the negative eigenvalue to the fluctuation 
along the real axis in Eq. (|38D . After performing the integrations we obtain the following 
overall normalisation factor (using ( [A.4| )) 



N{t,Y) 



1 to 
(27r)2cu 



(27r) 



3/2 _ 



buj 



a/ — A1A2A3 



-Anx'il) ^ 



(A.6) 



From (|A.6|) we can derive the subleading ~ Y"^ corrections. To this aim we expand, 
using (^3D, LJ and x'il) iii C as follows 

iU = cu,(t)(l+D^) + ... 
Xil = l-P) = -Px"il/2)-l^fx^'\l/2) = 



2^ 12 ^ ' 



x"(l/2) + ^(0'x(^Hl/2) 



+ 



(A.7) 



with D = Xm/C^Xm)- Inserting ( [A.7| ) into ( |A.6| ) gives 

5 



N{t,Y) 



+ 1dc')(i - -DC')(1 - -^)(i + -^C') 



1 



AttDusY 



1 + 



1 (4) 

1 Xm'Xm 



12 24 



(A.8) 



Appendix B: Calculation of and Y terms in the 
Airy approximation 

We show here how to obtain the quadratic ~ Y'^ and linear ~ Y corrections to the 
saddle point exponent uJs{t) in the case of the Airy diffusion model. We start from the 
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exponent in the expression (BT 



E 



^s{i)\Y\ 



(-^)/3 



+ 



Aty ^ 



2i J 



(B.l) 



and consider the expansion around the saddle point ^ defined by equation (|6^). We 
obtain 



E = ao + aiAe + ^a2{A^Y + ^a^iA^f + 



where we defined 



and 



(B.2) 



(B.3) 



(B.4) 



For simphcity we take = 0. The saddle point condition ai = results in equation 



([6^ ) which at small values of z has a solution (|63| ) ^ = 1 — . The term oq instead 
provides the leading exponent Ug. The second derivative evaluated at ^ reads then 



a2 



2uj,Y 



-3(3^-2) 1 
z e 



and the third one 



as 



3^ 1 



We then expand the exponential in the following form 



62^ 



(B.5) 



(B.6) 



(B.7) 



By integrating over the expression (B^) we arrive at the following result 

15 2^al 




— 1 + 



16(3!)2a3 



(B.8) 



The first term in (|B.8| ) i.e. the factor a/^ttTo^ comes from the second order fiuctu- 
ation and together with normalisation in ( |6T] ) gives the following overall normalisation 
factor to Gpert{y',t,to) and the term 



27r^D\z\/3\ ^ V 



1 ^ 



(B.9) 
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which checks with the first term ~ DQ"^ in Eq. ( [A.8| ) (due to the Airy approximation 
we will not reproduce the second term with higher order derivatives). 

The second term provides in turn the linear term in Y 



15 2^al 15 /3cj,r\ V z V 5 2 5 



16(3!)2a^ 16 V ^ y K^^sY ] 16 cu,F 32^ ' sAm , K j 

which is a second order shift to the saddle point value ujs{t). The (A^)^ term in ([B.7D 
can be safely neglected since it is of the subleading order (bas)'^. 
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